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We demonstrate quantum bath engineering for a superconducting artificial atom coupled to a 
microwave cavity. By tailoring the spectrum of microwave photon shot noise in the cavity, we create 
a dissipative environment that autonomously relaxes the atom to an arbitrarily specified coherent 
superposition of the ground and excited states. In the presence of background thermal excitations, 
this mechanism increases the state purity and effectively cools the dressed atom state to a low 
temperature. 



In practice, quantum systems are never completely iso- 
lated, but instead interact with degrees of freedom in the 
surrounding environment, eventually leading to decoher- 
ence of some states of the system. Precision measure- 
ment techniques such as nuclear magnetic resonance and 
interferometry, as well as envisioned quantum schemes 
for computation, simulation, and data encryption, rely 
on the ability to prepare and preserve delicate quan- 
tum superpositions and entanglement. The conventional 
route to long-lived quantum coherence involves minimiz- 
ing coupling to a dissipative bath. Paradoxically, it is 
possible to instead engineer specific couplings to a quan- 
tum environment that allow dissipation to actually pre- 
serve coherence 1-3 . In this letter, we demonstrate such 
quantum bath engineering for a superconducting arti- 
ficial atom coupled to a microwave frequency cavity. 
Cavity-assisted cooling of the atom is tailored to pro- 
duce any arbitrary superposition of ground and excited 
states on demand with high fidelity. 

The concept of our experiment is shown in Figure 1. 
A two- level atom is driven resonantly at frequency w q . 
In the frame rotating at the drive frequency, the two 
eigenstates of the system are |±) = (\g) ± |e))/v2, with 
eigenvalues ±1 of the a x Pauli operator. The energy 
splitting between the |+) and |— ) states is given by the 
Rabi frequency, Or. If M7r <C ^B^eff, where T e g is the 
effective temperature, neither state is thermodynamically 
preferred. However, by weakly coupling the atom to a 
cavity and introducing a drive detuned from the cavity 
resonance by A c = w,j - oj c , the photon shot noise of 
the cavity forms a quantum bath for the atom which can 
be engineered such that dissipation drives the atom to 
the |+) or |— ) state. Here, uj d (uj c ) is the drive (cavity 
resonance) frequency. For red detuned drive, (A c < 0), 
the cavity dissipation "cools" the atom to the |+) state. 

Cavity assisted cooling has been studied extensively 
in the context of atomic gases 4-7 , mechanical objects 8-11 
and spins 12 . Similarly, single atoms and qubits have been 
used to alter the dissipation environment of a resonator, 
leading to lasing 13,14 , cooling, and amplification 15,16 . We 
demonstrate that cavity cooling can be applied to the 
dressed states of a two-level atom, and the dissipation in- 
troduced by the drive may be engineered to relax the sys- 
tem towards any specified point on the Bloch sphere — a 



valuable resource in quantum information processing. 
This process is resonant and can produce large cooling 
rates. 




FIG. 1: Cavity cooling of a dressed state. A two-level atom is 
driven resonantly to form new eigenstates, |+) and |— ) in the 
rotating frame. Driving the cavity at A c = — fi^ resonantly 
enhances the anti-Stokes process relaxing the system to the 
|+) state. 

Our two-level atom is realized using the two lowest 
energy levels of a superconducting transmon qubit 17,18 
with w q /27r = 5.0258 GHz, dispersively coupled to the 
TEioi mode of a 3D superconducting cavity with fre- 
quency u! c /2ir — 6.826 GHz and linewidth k/2tv = 4.3 
MHz. The qubit induces a state-dependent frequency 
shift on the cavity of — x a z where x/2-7r = —0.66 MHz is 
the dispersive coupling strength. Similarly, the qubit fre- 
quency undergoes a light (or AC Stark) shift depending 
on the intracavity photon number, n, with mean value 
denoted w q = w q + 2\n- The Hamiltonian for the qubit 
in the frame rotating at w q is, 

H _ — j?,^ -x<Jaa z , (1) 

where a' (a) is the cavity photon creation (annihilation) 
operator. Transition rates between the |+) and |— ) 
states are determined from Redfield theory 19 , r± = 
[4x 2 Snn{T^R) + S yy (T^R) + S' zz {Tn R )}/A, where S yy ~ 
l/(Ti) and S' zz ~ 2/(T v ) are the power spectral den- 
sities of noise orthogonal to the x axis in the rotating 
frame, and S nn (ui) = n/s[(«;/2) 2 +(a;+A c ) 2 ] _:L is the spec- 
tral density of photon number fluctuations in the cavity 19 
that characterize the quantum bath. T\ and T v are the 
energy decay and pure dephasing times for the qubit, re- 
spectively. For A c = 0, S nn is symmetric in frequency 
and corresponds to an infinite temperature bath. In this 
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case S nn causes dephasing of the qubit and can equiv- 
alently be described in terms of a fluctuating AC stark 
shift, or in terms of measurement induced dephasing 19 ' 20 . 
When A c ^ 0, S nn is asymmetric and corresponds to a 
bath with finite positive (or negative) temperature and 
can be used to cool (or invert) the qubit state 19 . As il- 
lustrated in Figure la, cooling takes place via inelastic 
Raman scattering of pump photons. The most efficient 
cooling to the |+) state takes place for A c = — f2 R where 
the anti-Stokes photons are on resonance with the cavity. 
At this point the net cooling and heating rates are, 



r_ = 
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respectively, where T 2 = (l/2Ti + t/%)- 1 = 10.6 ^s 
is the (lab frame) dephasing rate. In equilibrium with 
the cavity bath, the final qubit polarization is given by 
detailed balance. 

The state of the qubit was probed by measuring the 
phase shift of a microwave tone reflecting off of the cavity 
at the cavity resonance frequency. The reflected signal 
was amplified by a lumped-element Josephson paramet- 
ric amplifier (LJPA) operating in phase sensitive mode 
which allowed high fidelity, multi-state, single shot read- 
out of the qubit state. Our qubit sample exhibited ex- 
cited (and higher excited) state population in excess of 
what would be expected from the nominal T — 20 mK 
environment. For our measurements, we used post se- 
lection to disregard the higher excited state populations 
which was as much as 12% of the qubit population. 



(a) 



ji/2 ji/2 
delay time 




0.000 0.002 0.004 0.006 0.008 0.010 

delay time (ms) 

n/2 




0.000 0.002 0.004 0.006 0.008 0.010 1.002 

cooling time (ms) 

FIG. 2: (a) Ramsey measurement using detuned pulses, (b) 
Cavity cooling to the |+) state with — A c /27r = Qr,/27T = 9 
MHz, and n = 3.6. After driving the system at uu'^ and uid for 
a variable time, a detuned n/2 pulse transfers the remaining 
coherence to the a z basis. High contrast, persistent "Ramsey" 
fringes indicate that the qubit has been cooled to the |+) 
state. 

To demonstrate effective quantum bath engineering, 
we compare a Ramsey measurement (Fig. 2(a)), to an 



experiment in which the qubit was cooled to the |+) state. 
The Ramsey measurement consisted of two 7r/2 pulses 
detuned by 2.8 MHz from the qubit frequency, followed 
by state readout in the a z basis and showed a typical 
T£ = 4.9 /Lts exponential decay of coherence. In Figure 
2(b), the cavity was driven to cool the qubit to the |+) 
state. After cooling the state for a variable period of time, 
the remaining coherence was measured by applying a tt/2 
pulse at a frequency detuned by 2.8 MHz from and 
measuring the amplitude of the resulting oscillations in 
the ground state population. We note that there is an 
initial build up of the coherence over a time scale of less 
than 1 fis given by = (r + + T_) _1 after which the 
system enters a steady state and the coherent oscillations 
persist indefinitely. 

In Figure 3 we present tomography of the qubit state 
after it has come to equilibrium with the dissipative en- 
vironment presented by the cavity. The qubit was driven 



at a variable detuning A c 



u r , where w r is the 



drive frequency, and at variable cavity drive power de- 
tuned from the cavity by A c /27r = —9 MHz. The am- 
plitude of the qubit drive was fixed to give f2ft/27r = +9 
MHz on resonance. In Figure 3(a) we display the tomog- 
raphy data for (a x ) and (a z ) (note that in the rotating 
frame (a y ) = 0). The dashed line indicates the depen- 
dence of ui q — Wq on the drive power, P^, indicating the 
detuning where the drive maintains resonance with the 
light-shifted qubit. Following this curve, we plot (tj x ) 
and (<t z ) as a function of cavity photon number in Fig- 
ure 3(b). When n = 0, T + = T_ = 1/(2T 2 ) and the 
system is completely incoherent. As the number of intra- 
cavity photons increases, coherence builds up along the 
x direction and saturates around n = 1. The purity of 
the cooled state is given by r_/(r + -fT_). The observed 
maximum state purity was 70% and was limited by our 
state readout fidelity (~ 90%) and population relaxation 
to the second excited state of the transmon. The latter 
reduced the measured state purity by up to an additional 
~ 20% depending on the time delay between the tomog- 
raphy pulses and the readout. The combination of these 
two effects led to a 70-80% reduction in the measured 
state purity. Taking into account these reductions, our 
measurements are close to the predicted value, plotted as 
a dashed line in Figure 3b. 

When |xV"| > K , the system is in a regime of strong 
coupling where higher-order and rapidly-rotating terms 
that have so far been neglected become significant. To 
explore this regime, we performed numerical simulations 
of the master equation. These results are shown in Fig- 
ure 3(b) as solid blue and gray lines and indicate that 
the maximum state purity is reduced compared to the 
predictions of Eq. (2) at high drive powers. Based on the 
simulation, we estimate that the actual state purity was 
94%, which corresponds to an effective temperature in 
the rotating frame of T e g — 150 /iK for the dressed state. 

To quantify the strength of the cavity damping, we 
plot the measured cooling rate in Figure 3(c), obtained 
by measuring the exponential timescale for buildup of 
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FIG. 3: State tomography, (a) Color plots show (a x ) and (a z ) as a function of cavity drive power (P<j = 101og(n) (dB) ) 
and qubit drive detuning for fixed cavity drive detuning, —A c /2tt — fl-n/2-K — 9 MHz. The dashed line indicates the optimal 
detuning A q = oj q — uj r = 0. (b,c) (o x ), (a z ) and the cooling rate are plotted versus n for to? = w q . The dashed lines indicates 
the prediction from Eq. (2) which has been scaled by our measurement fidelity of 80%. The solid blue, gray, and green lines 
indicate the results of the simulations for {cr x }, and the cooling rate respectively. Error bars in (c) represent the estimated 
error in the exponential fit. (inset in c) (a x ) vs. cooling time for n = 1.4. (d,e) Bloch sphere diagrams indicate that when the 
qubit drive is off resonance, the |±) states are tilted from the equator of the Bloch sphere, (f) Color plot shows inversion of 
(a x ) and (a z ) for +A c /2tt = fi R /27r = 9 MHz. 



ensemble population in the |+) state as the duration of 
the cooling pulse was increased (Fig. 3(c), inset). The 
measured rate is in quantitative agreement with Eq. (2), 
shown as a dashed line, as long as IxV^I < K - At higher 
photon numbers, the observed increase in coherence was 
not exponential. In this regime, the system is expected 
to exhibit damped oscillations between |— ) and |+) in 
analogy with vacuum Rabi oscillations (see the supple- 
mentary information). 

When the qubit drive is off-resonant, the engineered 
dissipation drives the qubit to different points on the 
Bloch sphere. As we illustrate in Figure 3(d), the off- 
resonant qubit drive creates an effective magnetic field 
at an angle 9 — arctan(f2^/A q ) with respect to the z 
axis, tilting the |±) states from the equator of the Bloch 
sphere. Here, A q = uj' q — u x is the detuning of the AC 
stark shifted qubit frequency from the qubit drive. In this 
case the cavity dissipation drives the system to the state 
obeying <Jg\+) = (+1)1+) where erg = sin 9a x + cos 9cr z . 
When the cavity drive is very weak (Fig. 3(e)) dissi- 
pation due to the finite qubit T\ = 10 ^s favors the 
|— ) state when A q > 0. In Figure 3(f) we display 



qubit state tomography when the cavity drive detuning, 
+A c /27r = 51r/27t = 9 MHz. In this case the cavity 
dissipation inverts the qubit to the |— ) state. 

By driving the qubit off-resonance and altering the cav- 
ity detuning to remain equal to the off-resonant Rabi fre- 
quency, J7r = yj f2p + Aq 2 , arbitrary superposition states 

of \g) and \e) can be prepared using the cavity dissipa- 
tion. As we show in the supplemental information, the 
heating and cooling rates are reduced by (^Ft/A q ) 2 since 
the qubit drive is no longer resonant. In Figure 4 we dis- 
play measurements that demonstrate cooling to arbitrary 
latitudes on the Bloch sphere. Figure 4 displays (ag) for 
9 = {90°, 43°, 10°} for A c /2tt = -15 MHz and variable 
drive power and detuning. For 9 = it/2 the cooling is 
along x as in Figure 3(a). As 9 is decreased (Fig. 4(b-c)), 
the optimum cooling occurs for A q = yA| — 

For a weak Rabi drive strongly detuned from w q , 9 ap- 
proaches zero, |+) ~ \g), and the cavity dissipation mech- 
anism crosses over to ordinary cavity-assisted Raman 
sideband cooling transitions similar to the transitions 
used to cool atoms 21-24 and superconducting qubits 25 us- 
ing a strong atomic transition to enhance emission at a 
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FIG. 4: Preparing arbitrary superposition states 
using cavity dissipation. Measurements of (ag) where 
o g = sm9a x + cos9a z for 9 = {90°, 43°, 10°}, (a-c), ver- 
sus drive power and qubit-drive detuning. The dashed lines 
indicate — oj q . The transparent gray lines indicate the 
qubit drive detuning that gives the most efficient cooling, 
(A 2 — fin) 1 / 2 . For 9 = 10°, two sideband transitions are vis- 
ible corresponding to Raman transitions that leave the qubit 
in the ground (lower inset) or excited state (upper inset). The 
ground state Raman process (lower inset) involves absorption 
of a cavity drive photon, stimulated emission at cjq + A c , and 
emission into the cavity at cj c . Two possible processes are 
labeled by gray and black arrows. 

specific frequency. Here, the cavity takes the place of the 
strong atomic transition. In Figure 4(c), with 9 = 10°, 
two "sideband" transitions emerge for Pd > dB. Level 
diagrams of the transitions are shown as insets in Fig- 
ure 4 indicating that simultaneously detuning the drives 



from the cavity and qubit allows selective optical pump- 
ing from the ground or excited states. The rates for these 
transitions can be calculated using Fermi's golden rule 
and agree with a calculation based on Redfield theory 
(see the supplementary information). 

As we previously noted, when A c = 0, photon num- 
ber fluctuations induce dephasing of the qubit in accor- 
dance with the theory of quantum measurement 19 . In 
this regime, cavity photons convey information about the 
qubit state encoded as a phase shift corresponding to an 
elastic scattering event. When the drive is detuned, the 
measurement is replaced by an inelastic scattering pro- 
cess in which the scattering of a photon into the cavity 
heralds a transition to the effective ground state of the 
system. 

In conclusion, we have demonstrated quantum bath en- 
gineering with a model two-level system. The technique 
allows arbitrary superposition states of the the system 
to be prepared simply with saturating pulses. The tech- 
nique does not require active feedback 26-29 and is not lim- 
ited by detection efficiency. State preparation fidelities in 
excess of 99.9% are in principle possible with currently 
achievable sample parameters (T2 = 150 /xs, k/2tt = 2 
MHz, tt R /2TT = 50 MHz, = -1 MHz). Future 

multi-qubit implementations could enable the prepara- 
tion of entangled many-body states suitable for quantum 
simulation and computation. 
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I. SAMPLE FABRICATION AND PARAMETERS 

The qubit sample was fabricated from aluminum deposited on high resistivity silicon in a two-step, double-angle 
evaporation process with an intervening oxidation step. The qubit consisted of a single Josephson junction with 
critical current Iq = 46 nA shunting two paddles that set the charging energy (Ec/h = 145 MHz) of the qubit and 
provide coupling to the cavity. The coupling rate of the qubit to the cavity, g — 70 MHz was set by placing the sample 
away from the cavity electric field antinodc. 

The cavity was machined from aluminum (alloy 6061) and had dimensions of 35.6 x 5.2 x 28.1 mm 3 . The cavity decay 
rate, k/2-k = 4.3 MHz, was set by adjusting the length of the center pin of a SMA coaxial connector protruding into 
the cavity volume. The cavity was also addressed by a weakly coupled port which was used for the qubit manipulation 
pulses. 

The dispersive coupling rate x/2f = — 0.55(.ll) Mhz was determined by measuring both the AC stark shift and 
measurement induced dephasing rate of a drive tone at the cavity resonance frequency as in our previous work 1 . 
We use x/27r = —0.66 MHz for comparison to theory and simulation because this value gave the best quantitative 
agreement. 

II. EXPERIMENT SETUP AND STATE MEASUREMENT 

Drive signals were combined at room temperature and were sent to the sample via heavily attenuated coaxial lines. 
The reflected signal from the strongly coupled cavity port was amplified by a near quantum limited lumped-element 
Josephson parametric amplifier (LJPA) operating in phase sensitive mode. The LJPA was separated from the cavity 
by 4 cryogenic circulators and isolated from the HEMT amplifier at 2.7 K by three more isolators. The LJPA was 
biased with a tone at 6.823 GHz and exhibited a gain of 19 dB with a 3 dB (instantaneous) bandwidth of 30 MHz. The 
phase of the LJPA pump was adjusted to amplify the quadrature of the reflected signal that contained information 
about the qubit states. 

State measurement was performed by pulsing the readout and LJPA pump tones on. After allowing 200 ns for 
transients to decay, a 200 ns section of data was integrated. After repeating the experiment many times, a histogram of 
the readout voltage revealed three well-separated distributions corresponding to the qubit ground, excited, and higher 
excited transmon states. From these measurements we determined that the equilibrium population of the qubit states 
was 77%, 14%, and 9% for the ground, excited and higher excited states respectively. During cooling, the population 
of the higher excited states increased to 12%, owing to the the increase in the average excited state population. In the 
experiment we use post-selection to remove the higher excited state population from our tomography measurements, 
but we note that coupling between the excited states results in a significant reduction of the measurement contrast. 
We measured a2/js timescale for equilibration between the excited states of the transmon, which results in a ~ 20% 
reduction in the state purity determined by tomography. 

We note that it is in principle possible to reduce the second excited state population using sideband cooling as 
demonstrated in Figure 3d simultaneously with cavity cooling to prepare pure states without the need to post select 
away population in higher qubit states. 

In a typical experimental sequence the qubit and cavity drives were suddenly turned on for a 50 /xs duration, followed 
by state tomography and readout. The experiment was repeated at 6 kHz. Each data point in the tomography plots 
was the result of 10 5 experimental sequences. 

III. DERIVATION OF THE SYSTEM HAMILTONIAN & MASTER EQUATION 

We begin our treatment with the Jaynes-Cummings Hamiltonian and two external drives on the cavity, one with 
amplitude ed and frequency Wd near the cavity resonance frequency w c , and the other with amplitude e r and frequency 
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w r near the frequency of the qubit cu q : 

H = uj c cJa - ^a z + g(a^a- + aa+) + [e d e- <Wd V + ^eT^a} + h.c] , (1) 

where a + is the operator exciting the qubit and cr_ de-exciting it. By including qubit and cavity decay terms we can 
write the master equation 2 : 

p=-i [H, p] + KV[a]p + Y -^V\a 7 \p + T_V[a_]p + T + V[a + ]p, (2) 

where T>[L]p = (2LpL^ — L^Lp — pL^L)/2, r v = ijr- is the phase relaxation rate and the population relaxation rate 
is given by jr = T_ + T + . 

By applying the dispersive shift, U = es( acr +~ a+<T -), where A — uj q — uj c , and moving to the rotating frames for 
both the qubit and cavity U c = e w+a " d *, U q — e~ l<7z ^ t , we can simplify the Hamiltonian. Taking terms only up to 
second order in £ we are left with (assuming e r and ed are real): 

H = -A c a f a - Aq ^ X a z - - x^aa z + e d (a t + a), (3) 

where, Or = — \ = = Wd - w c , A q = cj q — u T , and we have ignored terms rotating at w r — or higher 

frequencies. Note that the second-order correction also includes the well-known Purccll effect 3 which adds a term 
proportional to k to the qubit relaxation rate due to the coherent mixing of qubit and cavity excitations. 

We eliminate the drive term by displacing the field operator a = a + d. Note that this displacement acting on the 
cavity dissipator term introduces 

KD[a]p = nV[d]p + | [(a*d - ad!) , p] . (4) 
Because the last term has the form of a commutator, it is equivalent to a shift in the Hamiltonian of 

AH = y (a*d - a^) . (5) 



A convenient displacement choice will be: 



A c + in/2 



(6) 



With this choice we eliminate the drive and all of the terms linear in d in H + AH except for terms of the form da z . 
Our Hamiltonian can now be written as: 

H = -A c dU - Aq + 2 X(2n) a z - ^a x - X (a*d + aS + dU)a z , (7) 

where n = |a| 2 , and we have absorbed the Lamb shift into the definition of the qubit frequency [via the replacement 
[2n + 1) — > (2n)]. The Hamiltonian above is the one used in our simulations. 

To better understand this Hamiltonian and present a more intuitive understanding of the cooling process we choose 
the Rabi drive frequency such that A' q = A q + x(2n) = 0, meaning the qubit drive is exactly on resonance if there 
are n photons in the cavity. By performing a Hadamard rotation to interchange a z and <r x and writing the new a x as 
a + + O- , we arrive at an effective Hamiltonian (ignoring rapidly rotating terms as we choose the cavity drive (red) 
detuning to match the Rabi frequency, A c = —Or), 

H = -A c dU- ^a z - X (a*d<7+ +ad t cr_). (8) 

The new Hamiltonian is an effective Jaynes-Cummings model in which the effective qubit ground state is the lower 
energy eigenstate of a z [<r x in the original (rotating) spin frame before the Hadamard transformation]. With the 
Hamiltonian in this form, one understands that the cooling scheme is simply the dissipation of the dressed qubit to 
its effective ground state via photon emission into the cavity, just as in the ordinary (un-driven) Jaynes-Cummings 
model. [Here however the photon emission is actually Raman scattering of the cavity pump photons.] Using the fact 
that the peak value of the resonator density of states is p = ^ Fermi's Golden Rule yields for this particular resonant 
case, the following simple expression for the cooling rate 

r = 2^| X a| 2 p= ^-n, (9) 



which is valid in the limit of weak coupling (|x|\/n <§C ft). 
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IV. GENERAL HEATING AND COOLING RATES 



Before extending the above results to the more general case of arbitrary detunings for the cavity and qubit drives, 
it is convenient to reformulate the above derivation. Without making the Hadamard transformation and the cavity 
displacement transformation mentioned above, the Hamiltonian for the qubit alone can be written in the frame 
rotating at the Stark-shifted qubit frequency as: 



2x[a i a — n\ 

~2 



# q = — 7r<*x o ( 10 ) 



in which the dispersive coupling term can be thought of as adding photon shot noise that causes dephasing of the 
qubit 4,5 . Following the derivation of Fermi's Golden Rule in terms of noise spectral densities presented in 4 , the qubit 
excitation and de-excitation rates can be calculated from the spectral density of the noise perturbing the qubit in 
directions orthogonal to the x quantization direction. Spectral density at negative (positive) frequencies correspond 4 
to energy emitted (absorbed) by the bath: 

r± = \ {s„(TfiiO + s to (t^r)} , (ii) 

= \ {4x 2 S„ rl (TttR) + C(T^r) + S w (TfiR)} , (12) 

where S zz is the noise power spectral density related to the autocorrelation function of the noisy coefficient of a z in 
the Hamiltonian. It can be broken into two parts: one due to the photon number fluctuations, i.e., S nn and the other 
one due to all other processes S zz . T_ refers to the transition rate from the high-energy eigenstate of <r x with — 1 
eigenvalue to the low-energy state with +1 eigenvalue. This cooling transition is the dominant one for our assumption 
of a red-detuned cavity drive, i.e., A c = —Or,. 

Being cognizant 6 that S vv is evaluated in the rotating frame and with further approximations that: 

s w (TnR)^5 ww (o) = r 1) (13) 

C(T^r) ^ S zz (0) = 2r v , (14) 

5 "" M = ( K /2)2 + ( W + A c )2' (15) 

where T v is the pure-dephasing rate due to other sources, and we take S nn to be the result for an uncoupled driven 
cavity. Note that noise along the z direction is not affected by the transformation to the rotating frame. Using these 
terms we obtain 6 : 

,2 



For the cooling to be effective, we need: 



r+= (WT^T + T + T- (17) 



i.e., the shot-noise term in T_ should dominate so that the asymmetry in 5 , „„(±Or) strongly affects the transition 
rates. 

If the qubit is driven off-resonance by a detuning from its (Stark-shifted) frequency, the Hamiltonian takes the 
form: 

A^ n R 2 X Wa - n] 
# q = — - ~y°x g Gz ' ^ 



If the cavity detuning is set to the new Rabi frequency A c = ±y O^ + A^ 2 , it is possible to cool the qubit to an 
arbitrary position on the Bloch sphere using the same mechanism. The cooling and heating will occur between the 
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two states |±), which are the eigenstates of (A' q <7 z + Oro^) / J R + A^ 2 with eigenvalue ±1. 

The transition rates can be calculated in a manner identical to the on-resonancc case and arc found to be: 

I± = \ [S zz (r^l + A'^ sin 2 9 + S xx (r^l + A*) cos 2 9 + S yy (r^ 2 R + A^ 2 ) } , (20) 

~ X 2 S ml (T-0 R / sin0) sin 2 9 + sin 2 9 + ^ (l + cos 2 9) , (21) 

where we define tan# = Or/A^. 

With weak qubit drive at large detuning A' q >• Or, the qubit is no longer dressed and cooling is along the z 
direction. In this limit, the physics devolves into the ordinary sideband cooling process in which both the qubit and 
cavity drives are detuned by <~ A^. The rates are given by: 



(«/2)2 + 4A» \A'J 2 \4, 
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Note that the cooling is less effective than the resonant case, as we have an additional factor of f which suppresses 
the rate. 

V. RAMAN TRANSITION RATES 

To complete the picture, we provide a derivation of the transition rates for Raman sideband cooling to show that 
they indeed coincide with the Rabi cooling for A' q ^> Or. The Hamiltonian can be written as: 

H =H + V, (24) 

H =-^a z -A c dU-°^a x , (25) 

V = -x(ad t +a*d + d t d)a z . (26) 

Here we have displaced the cavity field a — a + d, with a = e<j/ (A c + in/2). The cooling/heating processes rely on the 
fact that the dispersive shift term couples the qubit and cavity. We thus treat the dispersive term as a perturbation 
and use Fermi's golden rule to calculate the transition rates. 
The qubit part of H can be further diagonalized and we get: 

|e}^|e) + ^| 5 >, (27) 
\9)^\9)-^r\e), (28) 



2A' 
q 



where we have used the fact that Or <C A q . Note that 



.,.l;r / ,.u| 2 = |(e ) 0|^ ) l>| 2 =X 2 n(^) (2!)) 

We also make use of the effective cavity density of states p(u) = — ^ Im ^_ A * +i!L . Applying Fermi's golden rule, we 
get the same cooling/heating rates as in the noise spectral density calculation, i.e., 

4 X 2 _/0 R ^ 2 



q 



r -=^ Af) ' ^ 



K-X 2 - /^R X 2 



r +=W2) 2+ 4AfUj ' (31) 
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FIG. 1: Comparison between theory and experiment (a,b) Color plots show (a x ) as a function of cavity drive power 
and qubit drive detuning for fixed cavity drive detuning, — A c = S7r = 2tt x 9 MHz for the experiment (a) and simulation (b). 
(c,d) Color plot of {a z ) for experiment (c) and simulation (d). Simulation results in (b) have been scaled by 0.8 to account for 
imperfections in the tomography. 




FIG. 2: Strong coupling time domain oscillations The figure shows the simulated time domain evolution of (a x ) (black) 
and (d'd) (grey). The total photon number in the cavity was taken to be n — 3.31 and k/2ty — 0.2 MHz. The oscillations 
between cavity photon and excited state qubit are clearly visible and the frequency of oscillation matches the theoretical 
prediction 2|x|V^ 3 the energy splitting in the effective Jaynes-Cummings model. 



VI. SIMULATIONS OF THE MASTER EQUATION 

A set of numerical simulations were performed to better understand our experimental results. All the simulations 
involved numerically solving the Master equation presented in Section III and especially the Hamiltonian in Eq. (7). 
The simulations were written in python and make significant use of the QuTip toolbox 7 . Figure 1 shows a comparison 
between the tomography results produced by the simulation and the experiment for (<j x ) in Figure la,b and (a z ) in 
Figure lc,d. The experimental results are seen to be consistent with the theoretical predictions of our model. 
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VII. STRONG COUPLING REGIME 



The theoretical calculations in the sections above are all based on the weak coupling assumption v^lxl ^ K > namely 
that the coupling constant in the effective Hamiltonian [e.g. in Eq. (8)] is much smaller than the decay rate. When this 
assumption starts to fail the calculated decay rates are no longer a good model for the behavior of the system, which 
instead of exponential decay to the effective ground state shows oscillations 8 in the time domain. Intuitively this can 
be described as a coherent oscillation of the system (following a sudden switch of the coupling) between two states. 
The first state has the qubit in the high-energy eigenstate of a x and no Raman photon in the cavity. The second has 
the qubit in the low-energy dressed state and one Raman photon in the cavity. The coupling is strong enough that 
the system coherently oscillates back and forth between these states multiple times before the cavity decay eventually 
brings the system to equilibrium. This regime can be reached experimentally with currently available parameters 
(strong drives and a low n cavity). Figure 2 shows a simulation of the system in the strong coupling regime in which 
this effect can be observed. One sees excellent agreement of the oscillation rate with the predicted effective 'vacuum 
Rabi frequency' 2|x|\/^. 
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